Form Factors and Decay Rates for Heavy A Semileptonic Decays 

from QCD Sum Rules 



R.S. Marques de Carvalho 1 , F.S. Navarra 1 , M. Nielsen 1 , E. Ferreira 2 and H.G. Dosch 3 
1 Institute) de Fisica, Universidade de Sao Paulo 
CP. 66318, 05315-970 Sao Paulo, SP, Brazil 

G\ 
ON , 

Q\ \ Institute de Fisica, Universidade Federal do Rio de Janeiro 

7-H ! 

Ih ' CP. 68528, Rio de Janeiro 21945-970, RJ, Brazil 

^Institut fiir Theoretische Physik, Universitdt Heidelberg 
Philosophenweg 16, D-6900 Heidelberg, Germany 

> 

cn ■ Abstract 
m 
O 

ON 
On 



43 

a: 



QCD sum rules for the determination of form factors of A& and A c semilep- 



Qh. tonic decays are investigated. With a form for the baryonic current appropri- 



ate for the limits of the heavy quark symmetries, the different tensor struc- 



tures occurring in the two- and three-point functions are separately studied, 

b : 

and in each case general relations are written for the form factors. Particular 
attention is given to the treatment of the kinematical region ascribed to the 
continuum. The t-dependence of the form factors and the decay rates are 
numerically evaluated and compared to experimental information. 

PACS Numbers : 12.38.Lg , 13.30.Ce , 14.65.Dw , 14.65.Fy 
1. Introduction 

In this paper we treat semileptonic decays of the heavy baryons A c and A&. The study of 
the matrix elements of weak decays of heavy hadrons are an important source of information 



on some elements of the Cabbibo-Kobayashi-Maskawa (CKM) matrix, but the determination 
of these fundamental quantities of the standard model require disentanglement from the 
effects of strong interactions occurring inside the hadrons. 

QCD sum rules |j] (for reviews see 0,0]) provide an important approach for calculating 
hadronic matrix elements and form factors for systems with both light and heavy quarks. 
This method deals with the nonperturbative aspects of QCD analytically using a limited 
input of phenomenological parameters and has been applied successfully to semileptonic 
decays of mesons. Unfortunately QCD sum rules are not expected to do as well for baryonic 
as for mesonic amplitudes. This difficulty is inherent to the method, as the basis of the 
sum-rule approach is an expansion in local operators. As a starting point the hadron is 
represented by a local interpolating field constructed from quarks. It is evident that for the 
case of baryons where at least three operators are taken at the same point this is a more 
drastic reduction than for the mesonic case with two operators. We therefore expect the 
contribution of higher resonances to be even more important for baryons than for mesons. 
Technically this reflects itself in the higher dimension of the interpolating fields and thus 
in a faster increase of the perturbative contribution, which makes the assumptions on the 
continuum contributions more decisive. 

In the limit where one of the quarks in the initial and final hadrons is infinitely heavy 
there are new flavor and spin symmetries, which have received much attention in the last 
years in the framework of the so-called Heavy Quark Effective Theory (HQET) @-^[ (for a 
review see |[L0|| ). QCD sum rules simplify considerably in the heavy quark limit but they 
can be applied also to evaluate 1/M corrections to the results obtained with infinite heavy 
quark masses. The method has had good success in the calculation of corrections to HQET 
(see jrD| and references therein). 

In the present paper, which is an extended and revised version of a previous letter [TT 



we evaluate the semileptonic decays of heavy A-baryons in the QCD sum rule approach as it 
was developed a few years ago for heavy meson decays []12|| . This approach treats full QCD 
but also reproduces the symmetries of HQET. In A c semileptonic decay, which is among the 



best investigated processes of this kind, the relevant non-Cabbibo suppressed CKM matrix 
element V cs is known, so that calculations of these decays may provide very good tests 
of the applied method. On the other hand, there are quite serious discrepancies between 
experiments and HQET in the ratios of the lifetimes of beautiful baryons to beautiful mesons 
(see e.g. [[TB|]), and hence it is of prime interest to investigate in full QCD all decay channels 
of the At baryon. 

We use the sum rule technique for three-point functions that was introduced |T^JT5|| 



for the study of the pion eletromagnetic form factor at intermediate Euclidean momentum 
transfer. We analyse our form factors in the physical region for positive values of the 
momentum transfer since in this case the cut in the t channel starts at t ~ m?Q (juq represents 
the heavy quark mass), and thus the Euclidean region stretches up to that threshold. 

This paper is organized as follows: in Sec. II we discuss the method and introduce the 
form factors to be evaluated. In Sec. Ill we collect and discuss the results for the two- and 
three-point functions sum rules in HQET and shortly discuss the results. In Sec. IV we 
present our results obtained from QCD sum rules for the form factors and decay rates for 
A c and A& decays, and finally in Sec. V we present our conclusions. Technical details about 
the method used are presented in Appendices A, B and C. 

2. Sum Rule Calculation 

We label generically the initial channel by / and the final channel by F. The weak 
transition current couples initial quark Qj (mass mj) to final quark Qp (mass mp), and 
we assume that the other two non-interacting quarks are massless (m u = = 0.). The 
initial and final baryon four- momenta are called pi and pp, with corresponding Mandelstam 
variables s — p\ and u = pp respectively. The current four-momentum is called q — Pi ~Pf, 
and t = q 2 . We also define the four- vector r = pj + pp. The initial and final baryons are 
called A j and A^, and their masses are M\ t and M\ F . 

In order to study the decay A/ — > A F + £ + v t using the QCD sum rule approach we 
consider the three-point function of the weak transition current from a initial to a final quark 
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= Qf1h{]- ~ 1§)Qi an d the interpolating fields of the initial and final baryons r}A T and 

Va f 

Mpf.Pi) =fj d'xd'y (0|T{r7 AF (x)J M (0)rJ A/ (y)}|0> e^V*" . (i) 

This three-point function of two spinors and a vector and axial vector can be expressed 
by a superposition of scalar form factors and the 24 vector valued 4x4 matrices of Table 1. 

The method of sum rules is based on the simultaneous evaluation of the expression 
of the correlator (|]) in a phenomenological representation obtained by inserting physical 
intermediate states (the usual spin | baryons and higher resonances), and in a theoretical 
representation obtained through the OPE expansion. This is detailed below. 



A. Phenomenological Side 

The phenomenological expression is continued by dispersion relations into the not-so- 
deep Euclidean region where pp < M\ F and p\ < M\ i and where it is believed that the 
theoretical representation is reliable. We introduce the couplings // and fp of the currents 
with the respective spin | hadronic states 

(0|t7a f |Af(pf)> = fFUpipp) , (2) 
(Ai(pi)\v Al \0) = fiMPi) i (3) 

and obtain the phenomenological representation of Eq. (|lj) 

n? \VF.Vl) = L v 2 _ M 2 (MPfWMPi)) 2_ M 2 

spins "F Ap Pi Aj 

+ higher resonances 

fF u f (pf) , , . u //5Tj(pj) 



+ higher resonances . (4) 

The general transition element (Ap(pp)\J fl \Ai(pj)) with a weak current between two spin 
7; states can be written in terms of the moments q = Pi — Pf and r = pi + pp, the spinors 
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uf{pf) and ui(pj) and the 24 structures made with 7 matrices (12 vector and 12 axial 
vector quantities) listed in Table I in the form 
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{Af(pf)\ MMPi)) = uf(Pf) E 



U < + ft < 



Ui(pi) 



(5) 



Table I 

Vector and axial vector structures constructed with four-vectors q, r and 
combinations of 7 matrices. Notation is from Bjorken-Drell, with a^ u = §(7^7 
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?M 75 
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7m 
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?"m #75 
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9m /75 
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^M /75 
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ia afl r a 




10 




ia Q/3 g a r /3 7 5 g M 


11 


i(T a pq a r l3 r tl 
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e p pcxpq p r p ^ a 



Since on-mass-shell A/ and baryons obey Dirac equations 

we may use the projector properties of the sums over the (on-mass-shell) spinors 

ui(pi)ui(pi) =j-pi + M Al , u F (p F )u F (p F ) = 7 ■ p F + M Af , (6) 

spins spins 

and the sum over 24 invariants in Eq. ([5]) can be reduced to only 6 independent terms. 
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The experimental information is usually represented |L6| by the invariant decomposition 
of the amplitude in the form 

(A J ,(p J ,)|J M |A J (p J )) = 

u(p F ) [ 7ll (F? + F A l5 ) + io» v q\FX + F A lb ) + q„{F 3 v + F A l5 )]u( Pl ) , (7) 

where the form factors F^' (i = 1,2,3) are functions of t = q 2 . 

We then obtain for the phenomenological representation of the correlator 

{JH>F + f' P 



n (phc„) ( 



Pf,PIj 



Pi - M\ F 
q,(F^ + F 3 A 75 



7„(*f + F A 7h ) + ia» v q v {F% + F A 75 ) 

Uih + fi) 



V\ - Ml 



+ higher resonances , (8) 



where we have defined 



f' F = f F M Ap , /; = fjM Kl . (9) 

The reason for the separate notations for / and /' is that they belong to terms with different 
spin content (with and without 7 matrix) and this will help in the classification of types 
of traces to be evaluated. Another reason is that we shall determine / and /' from the 
mass sum rules for Aj and A^, and there is no guarantee that we have exactly fj = f^M^, 
f' F = IfM Af . 

Referring to the spinor forms that appear in Eq. (|S]), we may pick up terms with or 
without 7 matrices in the numerators of the two fractions, and thus form four different 
kinds of products that we identify as 

F?> A (q 2 )fif F , F^ A {q 2 )M' F , F^)f t f F , F^tfWfr , (with i = 1,2,3) . (10) 

We project out a sum rule for each one of these 24 products by performing appropriate 
traces using Eq. @. Independent traces can be formed after multiplying Eq. (§) by the 12 
vector and 12 axial- vector structures of Table I 



= Tr [^rP' (phen) ] and u { = Tr [a;rP' (phen) ], % = 1, . . . 12 , (11) 



so that we obtain 2 x 12 independent relations that must be solved for the 2 x 12 unknowns 
quantities mentioned in Eq. ([TOD ( we m ust recall that before we must introduce q and r 
everywhere in the places of pi, pp)- The results are given in Appendix A. 

The constants //, fj, Jf, and f' F have to be calculated from two-point function sum 
rules, as usual. 



The theoretical counterpart is evaluated by performing the Wilson's operator product ex- 
pansion (OPE) of the operator in Eq. ([]]) and then taking the expectation value with respect 
to the physical vacuum. The term from the unit operator gives the usual perturbative con- 
tribution, while the vacuum expectation values of the other operators in the expansion give 
the nonperturbative corrections proportional to the condensates of the respective operators. 
Thus 



where the index i refers to the types and dimensions of the condensates. 

As usual in the sum rule method, we wish to evaluate the form factors that build the 
(Af\J,m\Ai) amplitude by matching the phenomenological representation Eq. (||) of the three 
point function with its theoretical counterpart from Eq. fli~2l) . The higher resonances are 
approximated by the perturbative contribution invoking quark hadron duality [0J. 

As it is well known from two-point sum rules for baryons [p!7]-[T9| , there is a continuum of 
choices for the interpolating currents. Of course the results should in principle be indepen- 
dent of the choice of the current (except for pathological cases which couple very weakly to 
the ground state), but the justification of the approximations depends on the choice made. 
From the point of view of heavy quark symmetries, the spin of the heavy A is carried by the 
heavy quark, with the light quarks being in a spin and isospin singlet state, namely 



B. Theoretical Side 




(12) 



r]A Q = e ABC (u A i 5 d B )Q' 



(13) 
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where u A and d B stand for the Dirac fields of light quarks of colors A and B, and Q c 
represents a heavy quark (b or c) of color C. This current couples strongly to the A states 
in the heavy quark limit, and with this choice the quark condensate and the mixed quark- 
gluon condensate do not contribute to Eq. (|i~2|). The gluon condensate does contribute, but 
experience with baryonic two-point functions and mesonic three point functions teaches us 
that it is of little influence. We are thus left with only the perturbative and the four quark 
condensate contributions. 

Once we have calculated explicitly the terms contributing to n^ theor - ) we can form the 
traces 

x, L = Tr [^rF' (theor) ] and u % = Tr [a;rP' (theor) ], i = 1, . . . 12 , (14) 
and thus project out the form factors given in Appendix A. 



C. Sum Rule 

To the 24 invariants in Eq. fl5[) correspond in the off-mass-shell case 24 invariant functions 
.Fj ' (s,u,t)fjfp (where by /) we mean either fj or f' T , and similarly for f F ) given in 
Appendix A. They obey double dispersion relations of the form 

ft^s, u, t) fp f p = i r ^ r du> , f + • • • , d5) 

7T Z Jm] Jm 2 F [S' — S){U— U) 

where — 4p^' A (s, u, t) equals the double discontinuity of the amplitude F^' A (s, u, t)f\ ' fp on 
the cuts rrij < s < oo, m 2 F < u < oo, which is given by the double discontinuity of the 
right hand side of the expressions in Appendix A. In Eq. (TTH) mi and mp are respectively 



the initial and final heavy quark masses and the dots represent subtraction polynomials 
in s and u, which will vanish under the double Borel transform \\L4\ , in a straightforward 
generalization of that used before |]. Applying the double Borel transform to Eq . fljjD and 
subtracting the continuum contribution we obtain 



f v,a 



71 Jrrij J nip 

xp^ A (s',u',t)e- s '/ M h- u '^ , (16) 



where Mf and M| are the Borel parameters and R CO nt(s', u', sq, uq) defines the continuum 
model with s and u$ being, respectively, the continuum thresholds for the Aj and A^ 
baryons, determined from the mass sum rules. 

We will consider two different region of integration for the continuum contribution. A 
rectangular region specified by a rectangular cutoff: 

Rcont(s, U, S , U ) = ©Oo - s')Q(u - u) , (17) 

and a triangular region determined by: 

R C ont{s',u',so,u ) =Q[u -—(s' - s ) -u) . (18) 

V 771/ / 

In the OPE side of the sum rules we can evaluate the double discontinuity of the traces 
Xi using Cutkosky's rules. The perturbative contribution to the double discontinuity of Xi 
in Eq.([L4]) is given by 

6 1 r{\Ts- m i) 2 r ^2 

VV [xi] (p er t) = 77; v7 — 1 I dm 2 m 2 / dK d\K\ d cos 9 K d(p K 

(27r) 3 8^A(s,«,t) J ° 1 



x5(K - K )5(\K\ 2 - \K\ 2 )5(cos6 K - cos0jf)8(l - coa9 K )Ni , (19) 



where 



and 



K = {s + m 2 -mf}/2y/s , \K\ 2 = X(s, m 2 , m?)/4s 



2p S0 K - u - m 2 + m 2 F 



where 



and 



2\Ps\\K\ 



p S0 = (s + U - t)/2y r s , \p" s \ 2 = X(s,U,t)/A£ 



X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . 

In Eq.(|TP|) Ni represents the traces of the Lorentz structure of the three-point function 
multiplied by each of the 12 independent vector structures t>f 
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Nt = Tr [(& -fC + mf)7 M (l - 75) (j>c ~ ft + miK] . 



(21) 



There are, of course, equivalent expressions for the double discontinuity of Ui (see Eq.([14])) 
where the Lorentz structure is multiplied by the axial vector quantities af . 
In order to estimate the four-quark condensate we use the factorization 



{d A u»u»d A ,} = — 2 6^6 Qa ,5 AA 5 BB (uu)(dd) 



(22) 



where the parameter k, which may have some value from 1 to 3, is introduced to account 
for deviations from the factorization hypothesis @. 

For the particular choice of current shown in Eq. fll3|) and using only the sum rules based 
on // and fp structures (for which the imaginary parts are positive definite) we obtain the 
general results 



F A (t) = -F?(t) ; F*(t) = F 3 v (t) = F A (t) = F A (t) = 



(23) 



with the t-dependence given by 



16 Jmj Jm% [\{s',u',t)f/ 2 
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&(u' - u min )&(u max - u') R cont (s' } u', s , uo) + -^K(qq) 2 e- m V M h- m V M l 



, (24) 



where 



P(s, u, t) = -{m^u + s[m A F + m 2 F (s — t — u) +tu] — Tn 2 [(s + t — u)u + m F (s — t + u)}} 2 , (25) 
and the quantities 



u 



2m 



Y ) m F (s + mj) + (s — rrij) 



mj-tT \l X(mj,mp,t) 



(26) 



define the physical region. 

As mentioned before, fi and fp are determined from the two point sum rules and are 
given by p| 
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p M i I{F) / M MHM F) 

h ^ ~ 207T 4 



5 m* r iU0) ds h(s) e-Wnxn 



,2 

■1(F) 



+ — (g 2 G 2 ) r (U0) (1 - x)(l + 5x 2 ) e - s/M «n^) + —n*K(qq) 2 e- m '(W) 
768 Jm 2 r/ „, 3 



/m 7(F) 



(27) 



with 



/ 8 1 \ 

/i(s) = (l-x 2 ) 1 -- + — )- 12 lnx, (28) 
V x ar/ 



where 



x — n?i(F^I s , 

while M MI / MF ) represent the Borel masses in the two-point functions of Aj and re- 
spectively. In Eq. (p7D we have included the gluon condensate, (g 2 G 2 ), contribution for 
completeness, but its contribution is negligible and does not influence meaningfully the 
actual calculation. 

Comparing Eqs. ( p4[ ) and (B7|) we can see that the exponentials multiplying Eq. ( p4|) 
disappear if we choose 

2M 2 MI{MF) = Mj {F) . (29) 

Indeed, this way of relating the Borel parameters in the two- and three-point functions is a 
crucial ingredient for the incorporation of the HQET symmetries, and leads to a considerable 
reduction of the sensitivity to input parameters, such as continuum thresholds Sq and uq, 



and to radiative corrections 20 



3. Leading HQ limit 

In the HQET limit for infinite quark masses, the expressions for the sum rules simplify 
considerably, as was first noted by Shuryak M . All form factors can be described in terms of 
a simple function, the Isgur-Wise (IW) function, £, normalized to 1 at maximum momentum 
transfer. We use the familiar ^/-variable, related to the square of momentum transfer by 

t = m 2 + m 2 F — 2mjm F y . (30) 

11 



For convenience we give the relation between the Dirac-type form factors used here and the 
velocity form factors appropriate for HQET. The latter are defined through 



u 



(pf) (kGY + vfcX + < Gl + l5 (j,G A + v^G A + < G A )) u( Pl 



(31) 



where v 1 and v F are the velocities of the initial and final hadron respectively, and the form 
factors are functions of y. 

Neglecting the matching conditions one has G\ = Gf = £(y)+o(£), and all other 
form factors are O ( 1 A . The relations with the form factors defined in Eq. (|7]) are given 



by [10,16 



F?(t) 

F 1 A (t) 
F 2 A (t) 



G\ + (m F + mi) 



1 G v + J_ G v 



2m 



2m j 



1 -GX - —G v 



2mi 



2m i 



3 > 



2mj 



1 G v + J_ G v 



2mp 



3 > 



-Gf - (m F - mi) ( J—G A + r— — G A 



\2m_f 



2mp 



1 G A + J_ Q A 



2mi 
1 

2m j 



2 



2m p 
1 

2mw 



3 > 



Go 



(32) 



The baryonic two-point functions were discussed both in full QCD and in the HQ-limit 
in ref. |19|] , three-point functions in ref. [21], and 1/m corrections have been considered in 
ref. Wll. 



In this section we collect and discuss the results for the two- and three-point sum rules. 
We stay on the firm grounds of Wilson's OPE and only use local condensates. The four-quark 
condensate, which is the most important non-perturbative input, is normally reduced to the 
square of the chiral symmetry breaking two-quark condensate through the factorization 
hypothesis. We allow for an uncertainty in the factorization hypothesis introducing a factor 
k that takes some value between 1 and 3 (see Eq. ([22]) ) . 
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Usually one introduces the variables 10 and /i related to the relativistic momentum square 
s and the variable M 2 which appears after Borel improvement by 

s = m 2 + 2mu ; M 2 = 2m/i , (33) 

where m represents the heavy quark mass (rrij or nip). Then the square of the coupling to 
the A current structure) is given by the sum rule derived from the two-point function 

h ( Js(/i) + ^ 2 s G 2 )h^) + yvrV(gg) 2 ) + O (1) , (34) 
with 

IJjj) = r° u n e-" /tM (ku . (35) 

J 

Here the cutoff uj c is the continuum threshold, i. e, above this value the resonances are taken 
into account by perturbation theory. One can easily see that Eq. (|34]) is obtained from 
Eq. fl27| ) when the limit m^p) —> oo is taken and the relations of form factors in Eq. ( |32| ) 
are used. 

In order to ensure Luke's theorem |23] for the three-point function one has to ensure that 
the two-point sum rule yields a correct difference between the baryon and quark masses. 
Using the reasonable values mj = 4.7 GeV and m c =1.35 GeV we obtain 0.94 GeV for the 
mass difference when the cutoff and Borel parameters take respectively the values u c — lA 
and n — 1 GeV or oj c = 1.3 and /x = 2 GeV. 



exp 



Ml 



m 



2m/i 
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FIGURES 
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FIG. 1. Quark-hadron duality in HQET. The physical region for the three-point function with 
heavy hadrons for a fixed value of y is the wedge between thesolid lines. According to ref. [25] 
duality can only be assumed after integration over the variable uj = uj\ — namely along the 
dotted lines. 



In order to determine the Isgur-Wise (IW) function, the choice of the continuum model is 
very important. As has been argued by Neubert |24| and shown in detail in a non-relativistic 
model by Blok and Shifman it is essential to use a triangular cutoff in order to take 
into account correctly the contributions of the higher resonances, since a rectangular cutoff 
yields non-analytic behavior of the IW function at y = 1 and, as has been shown at least 
in the model consideration in ref. p5| , quark hadron duality for the three-point holds only 
after the double spectral function is integrated along the dotted lines in Fig. 1. Performing 
this integration along the lines u = (oui + uj 2 )/2 = const, we are left with the sum rule for 
the three-point function, depending only on the Borel variable /i corresponding to u, and 
the IW function is given by 



/ 5 ID \ _1 



(36) 
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As discussed before, we have chosen the nonrelativistic Borel parameters for the three- 
point function to be twice the value of the two-point channel (here we have assumed the 
Borel parameters to be the same in the initial and final channels); this condition is essential 
for correct normalization. It can be shown that Eq. (p^ ) reduces to Eq. fl36|) when the limit 
m i{F) —> oo is taken. 
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FIG. 2. Isgur-Wise function at y=1.5 as function of the Borel parameter \i. 

1.0 




FIG. 3. IW function (solid line) and its derivative (dashed line) as function of y for [i = 1 GeV 
and u r = 1.4 GeV. 
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In Fig. 2 we show that the IW function at y = 1.5 as a function of the Borel parameter 
\x in the range from 1 to 2 GeV presents a reasonable stability. In Fig. 3 we display the IW 
function £(?/) and its slope £'(y) in the phenomenologically interesting range for y from 1 to 
1.5. The influence of the choice of the parameters can be read off from Table II, where we 
show the values of £(1.5) and of the slope £'(1). 



Table II 



Values of the IW function £ at y = 1.5 and of its slope £' 
at y = 1 for different values of the parameters. 





(GeV) 


(GeV) 


£(1.5) 


£'(1) 


standard set 


1 


1.4 


0.638 


-1.098 


standard set 


2 


1.3 


0.630 


-1.119 


standard set 


1 


1.6 


0.592 


-1.242 


no gluon condensate 


1 


1.4 


0.609 


-1.199 


no four-quark condensate 


1 


1.4 


0.560 


-1.333 


no four-quark condensate 


1 


1.15 


0.599 


-1.202 



We see that the strongest influences come from the choices of the cutoff u c and of the 
four-quark condensate, and also that the influence of the gluon condensate is small. 

We can represent the sum rule result for the IW function accurately through a pole fit 

^ = l + 1.136(y-l) • (37) 

It is important to remark that the slope of the IW function varies appreciably in the 
range 1 < y < 1.5 as can be seen in the dashed line in Fig. 3. Our result for the slope at 



maximum momentum transfer is about twice the value obtained in ref. 22 where a linear 
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fit is used. The slope at maximum momentum transfer is potentially important for model 
insensitive determinations of V c b [K§ from A b semileptonic decays. 
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FIG. 4. IW function as function of the momentum transfer t (see Eq. (|30|)) (solid line). The 
dashed line shows the same function with t and y related by the physical masses, as given in Eq. 



In order to compare with the results of the full QCD sum rules for the A& — > A c £ug decay 
we plot in Fig. 4 the IW function of Eq. ( pop as a function of t (solid line); the dashed line 
shows the IW function as a function of t when t and y are related with the physical hadron 
(and not the quark) masses through 



t = M 2 K + M 2 Ac - 2M Ab M Ac y . 



(38) 



4. Form Factors and Decay Rates in full QCD 

In this section we discuss the results for the structure fifp. The results based on the 
structures fjfp and fjf' F are given respectively in Appendices B and C. 

In order to determine the continuum thresholds Uq and Sq for the form factors and decay 
rates in full QCD we first analyse the mass sum rules. The expression for the two point sum 
rule in the ^-structure was given in Eq. (J27J) . The corresponding sum rule in the 1-structure 
is given by |L9 
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M luJ M MHMF) [ 



F) 



16vr 4 



8 mi ^ L 

O Jm I ( F) 



ds g(s)e 
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where x = mj^/s. 
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FIG. 5. Af, mass obtained from the mass sum rules, as function of the Borel mass squared for 
three different values of the continuum threshod so = (5.65 GeV + A s ) 2 : A s = 0.4 GeV (dotted 
line), A s = 0.5 GeV (solid line) and A s = 0.6 GeV (dashed line). 



In Fig. 5 we show M Ab as a function of the Borel mass Mub-, obtained by dividing 
Eq. (|D by Eq. @ for different values of 



s = (M Ab + A s 



(41) 



The parameter values used in all calculations are rrib = 4.6 GeV, m c = 1.4 GeV, M Ac = 
2.285 GeV, M Ab = 5.65 GeV and (qq) = -(230 MeV) 3 , and we have neglected the gluon 



condensate. As can be seen in this figure, there is a rather stable plateau for M^ b > 10 GeV 2 , 
for the three values of sq considered. A similar result is obtained for M\ c , with a stable 
plateau for M\ lc > 4 GeV 2 as shown in Fig. 6, where we plot M\ c versus the Borel mass 
Mmci for three different values of 

u = (M Ac + A u ) 2 . (42) 
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FIG. 6. A c mass obtained from the mass sum rules, as a function of the Borel mass squared for 
three different values of the continuum threshod uq = (2.285 GeV + A u ) 2 : A u = 0.6 GeV (dotted 
line), A„ = 0.7 GeV (solid line) and A u = 0.8 GeV (dashed line). The less stable result in the 
dot-dashed line is obtained for k = 2 (see Eq.(|22|)) and A u = 0.8 GeV. 

In Figs. 5 and 6 we have considered k = 1. A larger value of k decreases the value of 
Ma 6(c) and makes the curve less stable, as shown by the dashed-dotted line in Fig. 6, where 
we have used k = 2 and A u = 0.8 GeV. 

In ref. [|llj it was argued that the determination of the baryon mass based on the ratio of 
Eqs. ( |39| ) and ( |2~TD may be misleading especially for light baryons, since states with positive 
and negative parity contribute to it with opposite signs (we shall call this method 1). An 
alternative way to determine the baryon mass is based on the sum rule Eq. (^) and its 
derivative with respect to M^ 2 that yield 
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(43) 



with h(s) given by Eq. (j2"8|). Using the square root of the ratio between Eqs. ([43]) and 



(f27|) to determine F . (method 2) we have obtained curves very similar to the ones 
shown in Figs. 5 and 6, but with lower continuum thresholds. For instance, we cannot 
distinguish the solid lines in Figs. 5 and 6 from the new curves obtained with A s = 0.4 GeV 
and A u = 0.43 GeV respectively. Therefore, using method 2 we obtain similar stability as 
obtained with method 1, but lower values for the continuum thresholds, this effect being 
even bigger for A c . 



0.80 



0.60 



0.40 



0.20 



0.00 



0.0 



10.0 20.0 

M F 2 (GeV) 2 



30.0 



FIG. 7. Sum rule values for the decay amplitude at t=0 for the process Af, — > A c lv£ as 
function of the Borel mass Ml. The long-dashed line is the perturbative contribution, and the 
short-dashed line is the four-quark condensate for k = 1 (see Eq. (|2^)). The solid line represents 
the total contribution. 



We first discuss the process A;, — > A c £i>£. In the form factor sum rule result, Eq. 
introduce the relation between the two Borel masses 

Mj MX - nr, 
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In Fig. 7 we show the behavior of the contributions to the form factor F± at t=0 for 
the process A& — > A c £u£ for function of the Borel mass Mp , obtained by using 

a rectangular continuun cut-off, Eq. flI?D , in Eq. (|24|). We observe that the range in the 
Borel masses where the perturbative and non-perturbative contributions are in equilibrium 
is above 10 GeV 2 . Since the Borel masses in the two- and three-point functions are related 
by M^ T f MF -i = Mjr F J2 and given the relation in Eq. fl44|) , this means that Mp = 10 GeV 2 
is related with M^p ~ 5 GeV 2 and M 2 fJ ~ 15 GeV 2 which are reasonable values for the 
Borel parameter in the A c and A& mass sum rules, as can be seen in Figs. 5 and 6. To 
reproduce the experimental values of the A c and A& masses for these values of the Borel 
parameters we need A u ~ 0.65 GeV and A s ~ 0.45 GeV (method 1), which are the values 
used in Fig. 7. Using method 2 we obtain A u ~ 0.45 GeV and A s ~ 0.35 GeV which give 
a similar result for the form factor resulting from a smaller continuum contribution but a 
bigger four-quark condensate contribution. 

We have also calculated the t-dependence of this form factor in the range < t < 8 GeV 2 
(which covers the major part of the kinematically allowed region < t < 11.34 GeV 2 ) where 
we do not have difficulties caused by non-Landau singularities [T2|] . The t-dependence, 
represented with dots in Fig. 8, can be very well approximated by a pole fit, as shown by 
the solid line in Fig. 8 (method 1). The extrapolation of the fit to the maximal momentum 
transfer value, t max = {m^ — m^ c ) 2 , yields = —F^ = 0.97. In HQET this value is just 
the IW function at zero quark recoil (see Eqs. fl32])) and is predicted to be 1. 
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FIG. 8. Decay amplitude for the process — > A c £u£ as function of the squared momentum 
transferred to the leptons. Solid line: pole fit {t) = 8.46/(20.08 - t) to the sum rule results 



(dots) for as = 1 (see Eq. fl22|)) and A s = 0.45 GeV and A u = 0.65 GeV (see Eqs. (gg |42p) . Dashed 
line: the same for A s = A u = 0.5 GeV; F^(t) = 8.12/(22.27 - t). Dot-dashed line: the same for 
k = 2, A s = 0.5 GeV and A u = 0.8 GeV; F^ (t) = 9.84/(20.93 - t). 



In Fig. 8 we also plot the result obtained for F^(t) when using A s = A u = 0.5 GeV and 
M| = 15 GeV 2 as done in ref. fly] (dashed line). As can be seen in this figure, a very different 
value for F^(t mas[ ) is obtained in this case. The results obtained using A s = 0.35 GeV, 
A u = 0.45 GeV and M| = 10 GeV 2 (method 2) yield F^(t max ) = 0.84, with a t-dependence 
between the solid and dashed lines in Fig. 8, where we also show F^(t) calculated for k = 2 
(dot-dashed line). In this case the continuum thresholds that reproduce the experimental 
value of the A b and A c masses are A s = 0.5 GeV 2 and A u = 0.8 GeV 2 (see Eqs. ( fH"D and 
(H)) using method 1, and A s = 0.45 GeV 2 and A u = 0.55 GeV 2 (method 2). 

Another point already mentioned is the importance of the continuum model. The above 
results were obtained using a rectangular region of integration (see Eq. (0)). However, 
as pointed out before |2^,25|, to take correctly into account the contributions of higher 
resonances it is essential, in the HQ limit, to use a triangular region. Since in the present 



22 



problem the initial and final heavy quarks have different masses the triangular region cannot 
be the same symmetric triangle shown in Fig. 1, but it has to be modified by the ratio 
mi?/m/. Therefore, the limiting triangular region is defined by Eq. ([T8|). 
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FIG. 9. Same as Fig. 7 for the triangular region of the continuum model. 

In Fig. 9 we show the behavior of the contributions to the form factor F± at t — for 
k = 1 as function of the Borel mass Mp, obtained using the above defined triangular region 
of the continuum model and method 1 to determine the continuum thresholds. We observe 
that the perturbative contribution is now bigger than when using the rectangular region of 
integration and that the stability of the curve, as a function of the Borel parameter, is even 
better than in Fig. 7. 

The t-dependence of F^(t), extracted at Mp = 10 GeV 2 , is represented with dots in 
Fig. 10. Again it can be very well approximated by a pole fit as shown by the solid line. 
The dashed line shows the pole fit to the sum rule result for n = 2. The extrapolation of the 
fits to the maximal momentum transfer value, t max = (m^ b — m\ c ) 2 , yields FY = 1-05 for 
k — \ and Fi = 1.09 for k = 2 (method 1). These values have to be understood as upper 
limits to FY{t max ). 
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FIG. 10. Same as Fig. 8 for the triangular region of the continuum model. Solid line: pole fit 
FY(t) = 17.55/(27.94 - t) to the sum rule results (dots) for re = 1 (see Eqs. (§D . Dashed line: 
the same for k = 2; F^(t) = 20.57/(30.26 - t). 



Using the triangular region for the continuum model and method 2 to determine the 
continuum thresholds we get F^(t maiX ) = 1.01 for k = 1 and FY(t max ) = 0.99 for k = 2. 
Therefore, we conclude that the results are much less sensitive to the continuum thresholds 
when the triangular region of the continuum model is used. 

The slope of F{ as a function of y (see Eq. fl38|)) at maximum momentum transfer, that 
gives the slope of the Isgur-Wise function at zero recoil, £'(1), comes out to be very large 
for the rectangular region, £'(1) = —2.4 ± 0.4, but is more consistent with the results given 
in Table II for the triangular region £'(1) = — 1-45 ±0.15. The errors reflect variations of 
k from 1 to 2 and variations on the continuum thresholds obtained with methods 1 and 2. 
Also, comparing Figs. 10 and 8 with the dashed line in Fig. 4 we see that the results from 
HQET are very similar to the ones obtained with the triangular region in full QCD. 

In Table III we compare the HQET values for the form factors with 1/rriQ and short 
distance corrections given in ref. [IIJ with our results (see Eq. (|32|) ). The indicated errors 
account for both choices of models for the continuum integration. 
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Table III 

Form-factors for the semi-leptonic decay A& — * A c tvg at maximum momentum transfer for HQET 
in the limit of infinite quark masses, including I/thq corrections and short distance corrections 

pi and our results from QCD sum rules. 



form 



factors 



F? 



F v 
F v 



F A 
F A 

^3 



UlQ — > OO 



1 



-1 





with 
corrections 



1.03 
0.10 
0.00 
-0.97 
0.00 
0.12 



our results 



1.03 ±0.06 



-1.03 ±0.06 





With the pole parametrization of the form factors we can evaluate the decay rate of the 
semileptonic decay Aj, — > A c tvi. Using V<$> = 0.04, we obtain for the decay width 

T(A 6 -> Aj-Pt) = (3.3 ± 1.1) x 10~ 14 GeV , (45) 

where the errors reflect variations of k from 1 to 2, the choices of different regions of the 
continuum contribution and different continuum thresholds. This value is in agreement with 
other predictions [16,27-p9[| , which are in the range (3.5 to 6.17) xl0~ 14 GeV, and is also 



in agreement with the experimental upper limit [3(| given by 



r(A b -> A+£-v e + anything) = (4.4 ± 1.8) x 10" 14 GeV . (46) 
There is only one calculation of the semileptonic decay A& — > A c £u£ from lattice QCD PT 



which gives, however, only a lower limit to the decay rate since the method cannot access 
all the velocity transfer range. 
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FIG. 11. A mass obtained from the mass sum rules, as a function of the Borel mass squared 
for different values of the continuum threshod no = (1.115 GeV + A M ) 2 . The dotted line shows 
the result obtained from Eqs. ( pTp and fl39| ) and A u = 0.5 GeV. The other lines show the results 
obtained from Eqs. (|7|) and (||) and A u = 0.45 GeV (dashed line), A u = 0.5 GeV (solid line) and 
A u = 0.6 GeV (dot-dashed line). 

We now turn to the analysis of the A c — > Kivn semileptonic decay. We first analyse the 
A mass sum rule that can be obtained from Eqs. (|27|), (|39|) and (|43|) considered only up to 
order 0(m 2 s ). In Fig. 11 we show M\ as a function of the Borel mass Mms, obtained by 
dividing Eq. fl39|) by Eq. (|27| ) (method 1) for A M = 0.5 GeV (dotted line), where now 

u = (M A + A u ) 2 , (47) 

and m s = 0.17 GeV, Ma = 1.115 GeV. As can be seen from the dotted line in this figure, 
the A mass obtained by method 1 is not stable as a function of the Borel mass. In Fig. 11 
we also show M\ as a function of the Borel mass, obtained by dividing Eq. (^) by Eq. ( |27|) 
(method 2) for different values of A u . With this procedure we obtain a very stable plateau 
for M 2 Ms > 1.5 GeV 2 . 

Considering the masses of the baryons in the initial and final states the relation between 
the Borel masses in the three-point function is now given by 
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FIG. 12. The sum rule values for the decay amplitude at t=0 for the process A c — ► A£u£ as 
function of the Borel mass Mp for k = 1. The long-dashed line shows the perturbative contribution, 
and the short-dashed line is the four-quark condensate for k = 1 (see Eq. ([22])). The solid line 
shows the total contribution. 

In Fig. 12 we show the behavior of the contributions to the form factor F± at t=0 for 
the process A c — > A£i/£ for k = 1 as function of the Borel mass Mp . We observe that the 
region in the Borel mass where the perturbative and non-perturbative contributions are in 
equilibrium is above 4 GeV 2 . Since the Borel masses in the two- and three-point functions 
are related by M^j^ps = Mf, F J2, and given Eq. Q48|) , we have that Mp = 4 GeV 2 is 
related with M 2 MF ~ 2 GeV 2 and M\ n -5.5 GeV 2 . These are reasonable values for the 
Borel parameter in the A c and A mass sum rules, as can be seen in Figs. 6 and 11. To 
reproduce the experimental values of the A c and A masses for these values of the Borel 
parameters we need A n ~ 0.55 GeV and A s ~ 0.62 GeV (method 1) , which are the values 
that define the rectangular region of integration used in Fig. 12. Comparing Figs. 7 and 
12, we see that the importance of the four-quark condensate is much larger in A c than in A^ 
decay. Using method 2 to fix the continuum threshold for A c we get A s ~ 0.42 GeV and, in 
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this case, the importance of the four-quark condensate is even bigger, yielding a form factor 
with almost the same Borel mass dependence as obtained with nethod 1 but about 1.2 times 
bigger. 
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FIG. 13. The decay amplitude FY for the process A c — > A£u£ as function of the squared 
momentum t transferred to the leptons. Solid line gives the pole fit FY (t) = 2.52/(5.09 - t) for 
K = 1 and the rectangular region of the continuum contribution. The dashed line represents the 
same for k = 2, with FY (t) = 3.67/(7.07 — t). The dot-dashed line shows the same for k = 1 and 
the triangular region of the continuum contribution, with FY (t) = 4.25/(7.09 — t). 

The t-dependence of the form factor FY in the range < t < 0.6 GeV 2 where we do not 
find difficulties with non-Landau singularities |12[ is shown by the dots in Fig. 13 for two 
different choices of the four-quark condensate (k = 1, 2 ). The pole fits (solid line for k — 1 
and dashed line for k = 2) are also shown in this figure. For k = 2 the continuum thresholds 
that reproduce the experimental masses for M\ ls ~ 2 GeV 2 and M\ lc ~ 5.5 GeV 2 are 
respectively A u ~ 0.65 GeV and A s ~ 0.8 GeV (method 1) and A s ~ 0.52 GeV (method 
2). These differences in the continuum thresholds lead, in this case of the A c decay, to very 
different behavior of the form factors as a function of t. This change is not observed in the 
A b decay. 
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If we use the triangular region instead of the rectangular region for the continuum we 
obtain larger values for the form factor, as shown by the dot-dashed line in Fig. 13 for k — 1 
and A s ~ 0.62 GeV (method 1). In this case we can see that the pole file is not as perfect 
as in the other cases. For n = 2 and the triangular region the pole fit to the sum rule results 
is given by F?(t) = 5.34/(8.85 - t) (method 1). 

Using the form factors obtained with method 1 for the determination of the continuum 
threshold for A c and V cs = 0.975, we get for the decay width 

T(A+ -> A + e + + v e ) = (8.7 ± 1.2) x 10~ 14 GeV. (49) 

The overwhelming sources of the theoretical error are the uncertainties in the value of the 
four-quark condensate and the choice of the model for the continuum. Variations of the quark 
masses within reasonable limits such as 0.13 < m s < 0.17 GeV and 1.25 < m c < 1.45 GeV 
have negligible effects compared to the errors introduced by the variation of the four-quark 
condensate and the model of continuum. Within the errors the above result agrees with the 
reported experimental value [[!(J 

r cxp (A+ -> A + e+ + v e ) = (6.3 ± 1.9) x 10~ 14 GeV . (50) 

The form factors obtained with method 2 for the determination of the continuum thresh- 
old for A c yield a very large decay width ((13.6 ± 2.3) x 10 -14 GeV) but we do not think 
we can trust this result since, in this case, the sum rule is dominated by the four-quark 
condensate. 

Several theoretical attempts have been made to describe the form factors in the A+ — > 



A + e + + v e semileptonic decay, employing flavor symmetry and/or quark models f32|-|3~5" 
Their results for the decay width are in the range (4.7 - 10) xl0~ 14 GeV and are thus 
in agreement with our result in Eq. (fOJ. In spite of this good agreement in the decay 
width, Eq. ( pB"D leads to an asymmetry parameter a = — 1 while the observed value |30| is 
a = -0.82 ±0.10. 

We can also make predictions about the decay width of the semileptonic decay A;, — > piue. 
This decay was recently studied in the partial HQET framework |36| with result r(A& — > 
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p£ve) = (1.43 ± 0.07) x 1 11 1 V^fo | 2 GeV. Since the proton is not a heavy baryon, it is very 
important to compare this result obtained in the framework of partial HQET with our full 
QCD result though it should be remembered that there are serious objections against the 



use of conventional sum rules for heavy light transitions [[37] . 

The procedure is exactly the same as described above for the A c semileptonic decay. The 
results for the proton mass and form factor as functions of the Borel mass are very similar 
to those presented in Figs. 11 and 12, being only a little less stable. The region in the Borel 
mass where the perturbative and non-perturbative contributions to the form factor are in 
equilibrium is above 3 GeV 2 . Therefore, we evaluate the sum rules at Mp = 3 GeV 2 which 
give M\ 4p ~ 1.5 GeV 2 and M\ lh ~ 18 GeV 2 for the value of the Borel parameters in the 
mass sum rules of the proton and A& respectively. These seem to be rather large values for 
the Borel parameters in the two-point function, but we found that there is no stability in 
the form factor for smaller values of these parameters. 

To reproduce the experimental values of the proton and A& masses for these values of the 
Borel parameters we need A u ~ 0.6 GeV and A s ~ 0.4 GeV for k = 1 and A u ~ 0.75 GeV 
and A s ~ 0.5 GeV for k — 2, where A u and A s were determined by method 2 and 1 
respectively. Using also either the rectangular or the triangular region for the continuum 
contribution we obtain 

T(A b p + r + v t ) = (1.7 ± 0.7) x 10- u |y u6 | 2 GeV , (51) 



in agreement with the result of the HQET calculation ||36|| . 
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FIG. 14. Normalized charged lepton energy spectrum for the processes A c — > klvi (dot-dashed 
line), A& — > A c £i/£ (solid line) and — ► (dashed line). 

The knowledge of the form factors also allows us to calculate the charged lepton energy 
spectrum. In Fig. 14 we plot {l/Y)dY / dEi versus the energy of the charged lepton, Ei, for 
the three studied decays. The spectra have exactly the same shapes for k = 1 or 2 and for 
either the rectangular or triangular region in the continuum model. 

5. Summary and Conclusions 

We have developed in this paper the complete kinematical formalism to obtain the phys- 
ical form factors of semileptonic hyperon decays from the off-shell three-point functions and 
QCD sum rules. We have obtained the t -dependence of the form factors directly from the 
sum rules in the physical region for positive values of the momentum transfer since, in this 
case, the cut in the t channel starts at t ~ uiq (mg being the heavy quark mass), and 
thus the Euclidean region stretches up to that threshold. We have applied this formalism 
to the A+ — > A + e + + u e , A b — > A c + 1~ + Vi and A b — > p + 1~ + Vi semileptonic decays. 
It turns out that the sum rule for A+ — > A + e + + v e depends much more strongly on the 
values of parameters used in the sum rules than what has been found in the case of the 
analogous mesonic decay D — > + leptons . In particular, the dependence on the choice 
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of the model for the continuum makes detailed predictions impossible. This is not surprising 
since: a) the higher dimension of the baryon currents implies a faster increase of the dis- 
continuities along the cuts and thus a stronger dependence on the assumed onset of duality 
(continuum threshold); and b) spin leads to further complications. In our projection we 
have only singled out the ground state with spin 1/2 and put all remaining contributions in 
the pertubatively treated continuum. For higher spin states, which are also present in the 
three-point function, the projection is much more involved and such a simple treatment may 
be less justifiable than in the scalar and vector meson cases. This difficulty reflects itself 
in a certain inconsistency of the sum rules for the same form factor obtained from different 
gamma- matrix structures (see Sec. IV and Appendices B and C). We have therefore studied 
in more detail the fjfp structure (see Eq. (|24])). Another special feature of the A c semilep- 
tonic decay is the importance of the four-quark condensate which, for some values of the 
continuum thresholds, gives much larger contribution than the perturbative part. According 
to the original S VZ philosophy |I| , this brings doubt to the use of the sum rule approach in 
this case, and, furthermore, we must remark that the value of the four-quark condensate is 
only poorly known. So it is natural that the present analysis of the A c decay leads to very 
large errors, and the same is true, to an even larger extent, for A& — > p£~ve- 

The situation for — > h c i~vt is much more favorable. Here the symmetries valid in 
the limit of the masses going to infinity make themselves remarkable also in our treatment 
of full QCD. The result of the present more complete investigation modifies our earlier as- 
sertion |TT| about large 1/M 2 corrections for this decay. This is mainly due to the choice 



of the appropriate continuum model. It has been strongly argued by Blok and Shifman [25 
that only a triangular cutoff (see Sec. IV) is adequate for HQET, since a rectangular cutoff 
does not project out the p-waves properly. This argument becomes much more important 
for fermions than for mesons, since the change from a rectangular to a triangular cutoff, 
which had no strong effect on the decay rate of the mesons, is very important for baryons. 
So we conclude that a treatment of the — > KJrvn through HQET sum rules is satisfac- 
tory. Nevertheless, we must remark that the influence of the continumm threshold on the 
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phenomenologically important slope of the Isgur-Wise function is not negligible also in the 
HQET limit. 

We have calculated the semileptonic decay form factors using sum rules with three dif- 
ferent tensor structures. The most robust result that we have obtained for the form factors 
is that Ff(t) = —Fi(t) for all decays and structures considered. Different structures lead 
to different relations between F 2 y and F^- The only possible way to conciliate the three 
results would be to have F^ = F^ = 0, which however is a result valid only for the fifp 
structure. In the case of the A b decay, the numerical values obtained for F 2 V ^ in the f'jfp 
and fif' F structures are very small and, therefore, consistent with zero. However, in the case 
of the A c decay, the numerical values obtained for F 2 y ^ in the f' 1 f F and fjf'p cases are not 
consistent with zero, and this leads to huge variations in the decay rates obtained in the 
different cases. 

For the process A5 — > A c £vg our results for the decay width are in agreement with 
the experimental upper limit, and the extrapolation of the form factors to the maximal 
momentum transfer agrees very well with the HQET prediction. 

We have also calculated the charged lepton energy spectrum and it is interesting to 
remark that the form of the spectrum is not sensitive to the structure used, neither to the 
model of continuum, neither to the parametrization of the four-quark condensate. 
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Appendix 1 . Projecting out the Form Factors 

Defining 

^ = (uS-Ml) ' (Wlth i = W) ' (52) 
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and using 

A(s, t, u) = s 2 + u 2 + t 2 - 2su - 2st - 2ut , 
the 24 unknows denned in Eq. (|TUp are given in term of the traces by 



1 



F l (s,u,t)fif F = - 

A[X(s,t,u)] 2 

t\(s, u, t) x 3 + (2s 2 + 2st - t 2 - Asu + 2tu + 2u 2 ) x 4 
—3t (s — u) x 5 — 3t (s — u) xq + 3t 2 x 7 

F?(s,u,t)fifp = 1 • 
8[X{s,t,u)\ 2 

2(s — u)X(s, t, u) x\ — 2t\(s, t, u) x 2 — tX(s, t, u) x 8 + tX(s, t, u) x 9 
2s 2 — st — t 2 — Asu + htu + 2u 2 ) x w + 3 1 (s — t — u) x n 



FV( s ,u,t)f'J F 



8[X(s,t,u)\* 

2(s — u)X(s, t, u) xi — 2tX(s, t, u) x 2 — tX(s, t, u) x 8 — tX(s, t, u) x$ 

+ (2s 2 + hst-t 2 -Asu-tu + 2m 2 ) x w - 3 1 (s + t - u) x n 

F^(s,u,t)fjf' F = 1 • 
8[X{s,t,u)\ 2 

-t(s -t + u)X(s, t,u) x 3 - (s -t + u) {2s 2 + 2st - t 2 - Asu + 2tu + 2m 2 ) x a 
+3t (s - t + u) (s - u) x 5 + (s-u) (2 s 2 - st - t 2 - Asu - tu + 2 u 2 
—t (2 s 2 — st — t 2 — Asu — tu + 2 u 2 ^j x-j + itX(s, t, u) xu 
F 2 v (s,u,t)fjf F -- ' 



A[X(s,t,u)] 2 
X(s, u, t) x$ — 3 (s — u) Xiq + 3t xu 

F^(s,u,t)fif' F = - 1 • 
8[X{s,t,u)} 2 

(s — t — u)X(s, t, u) £3 + (As 2 — 5st + t 2 + Asu + tu — 8u~ ) x 4 

— (s 2 + st — 2t 2 — 2 su + 7 tu + u 2 ^j x 5 — 3 (s — t — u) (s — u) x 6 

+3 1 (s — t — u) x-j — iX(s, t, u) x\2 

1 



F^(s,u,t)f'J F 



l[X(s,t,u)f 
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— (s + t — u)\(s, t, u) rr 3 — (8s 2 — st — t 2 — 4sm + 5tu — 4m 2 ) X4 
+ (s 2 + 7 s t - 2 1 2 - 2 s u + t u + m 2 ) x 5 + 3 (s - u) (s + £ - u) x 6 
—3t(s + t — u) x 7 — ?A(s, t, u) X12 



F 2 y ( s ,M)/;/; = 



S[X(s,t,u)] 2 

— 2 (s — u) A(s, t, m) Xi + 2tA(s, t, m) x 2 + (s + m) A(s, t, m) 

(s — u) A(s, t, m) Xg + 3 (s — u) (s — t + m) X10 
(s 2 + s t - 2 1 2 - 2 S M + t U + M 2 ) X U 

1 



F 3 ( s , M ,t)M F = 4[A(s?t?M)]2 

A(s, m, t) x 9 — 3 (2 s — t + 2 u) x 10 + 3 (s — u) xu 



8[A(M,m)] 2 

(s — t + 3m)A(s, t, m) a: 3 + 3(2s — t + 2u)(s — t + 3m) x± 
-3 (s-u) (s - t + 3m) rc 5 - (5 s 2 -7 st + 2t 2 + 2su-tu- 7m 2 ) x 6 
+ (2 s 2 — st — t 2 — 4sm + 5£m + 2 m 2 ) x 7 — i A(s, £, m) X12 



8[A(M,u)] 2 

(3s — t + m)A(s, £, m) 0:3 — 3(3s — t + u)(2s — t + 2u) £4 
+3 (3s-t + u) (s-u) x 5 + (7 s 2 + st - 2t 2 - 2sm + 7tu - 5m 2 ) x 6 
— (2 s 2 + 5 s t — t 2 — 4sm — tu + 2m 2 ) x 7 + i A(s, t, m) x 12 

F?(8,U,t)f I f F = 



.[A(s,t,M)] 2 

2(2s - £ + 2m)A(s, t, u) xi + 2(s - m)A(s, t, u) x 2 
+ (s — m)A(s, t, m) x$ — (s + m)A(s, i, m) Xg 

+ (5 s 2 - 7 s t + 2t 2 + 14sm - 7tM + 5m 2 ) x 10 - 3 (s - t + u) (s - u) x n 
and for the axial form factors we obtain 
F^u^fjfp 



4[\(s,t,u)]< 
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t\(s, t, u) m 3 + (2s 2 + 2st - t 2 - Asu + 2tu + 2m 2 ) m 4 
—3t (s — u) u 5 — 3t (s — u) uq + 3t 2 u-j 
F^s^t)fjf F 



8[X(s,t,u)] 2 

— 2(s — u)X(s, t, u) U\ + 2tX(s, t, u) U2 + t\(s, t, u) u$ — t\(s, t, u) Mg 
+ (2s 2 - st-t 2 - Asu + 5tu + 2m 2 ) u w - 3 1 (s - t - u) u u 

Ff(s,u,t)f I f F = \ • 
8[X{s,t,u)\ 2 

2 (s — u) X(s, t, u) ui — 2tX(s, t, u) u 2 — tX(s, t, u) u 8 — tX(s, t, u) u 9 
+ (2s 2 + 5st - t 2 - Asu -tu + 2m 2 ) u 10 -3t (s + t-u) u n 



4X(s,t,u)\* 

t(s-t + u)X(s, t,u) u 3 + (s -t + u)(2s 2 + 2st - t 2 - Asu + 2tu + 2u 2 ) m 4 
—3t (s — u) (s — t + u) — (s — u) (2 s 2 — st — t 2 — Asu — tu + 2 m 2 ) m 6 
+t (2s 2 — st — t 2 — Asu — tu + 2m 2 ) u 7 — i tX(s, t, u) u 12 



F 2 A (s,u,t)f T f F 



1 



A[X(s,t,u)\* 
X(s, t, u) u 8 + 3 (s — u) u w — 3t uu 

Ff(s,u,t)f I f F = - \ • 
8[X(s,t,u)] 2 

(s - t- u)X(s, t, u) u 3 + (As 2 - 5st + t 2 + Asu + tu- 8m 2 ) m 4 
— (s 2 + st — 2 1 2 — 2 su + 7 tu + m 2 ) M5 — 3 (s — u) (s — t — u) u§ 
+3t (s — t — u) u 7 — iX(s, u, t) M12 
F 2 A (s,u,t)f'J F = 



^[X(s,t,u)} 2 

s + t — u)X(s, t, u) m 3 + (8s 2 — st — t 2 — Asu + htu — 4m 2 ) m 4 
— (s 2 + 7 s t — 2t 2 — 2su + tu + m 2 ) m 5 — 3 (s — u) (s + t — u) m 6 

+3 t (s + t — U) U 7 + iX(s, M, t) U\2 



8[X(s,t,u) 
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— 2 (s — u) X(s, U, t) U\ + 2tX(s, t, u) U2 + (s + u) X(s, t, u) u$ 

— (s — u) X(s, t,u) Ug + 3 (s — u) (s — t + u) U\Q 

- (s 2 + st- 2t 2 - 2su + tu + m 2 ) u n 

F£(s,u,t)f I f F = 1 • 
4[X{s,t,u)} 2 

— X(s, t,u) Ug + 3 (2 s — £ + 2 w) m 10 — 3 (s — w) uu 

^ M) ^ = 8[A(M, M )f - 

(s — t + 3u)X(s, t, u) U3 + 3(2s — t + 2w)(s — t + 3w) m 4 

— 3 (s — t + 3 m) (s — w) m 5 — (5 s 2 — 7st + 2t 2 + 2sw — tu— 7 w 2 ) m 6 

+ ^2s 2 — st — t 2 — Asu + 5iw + 2u 2 ^j u 7 — iX(s, t, u) m 12 



4X(s,t,u)} 2 

(3s — t + m)A(s, t, u) u 3 + 3(3s — t + w)(2s — t + 2m) u 4 
— 3 (3 s — t + u) (s — u) u 5 — (j s 2 + st — 2t 2 — 2 su + 7 tu — 5 m 2 ) u 6 
+ ^2s 2 + 5st — t 2 — 4su — + 2 ti 2 ^ M7 — iX(s, t, u) uyi 



3[A(s,t,u)] 2 

— 2(2s — t + 2u)X(s, t, u) u\ + 2(s — u)X(s, t, u) + (s — u)X(s, t, u) u$ 
(s + u)X(s,t,u) u 9 + (5s 2 -7st + 2t 2 + Usu-7tu + 5u 2 ^j u w 
3 (s — t + u) (s — u) u n . 



(72) 



(73) 



(74) 



(75) 



(76) 



Appendix 2 . Form Factors in the f'jfp Structure 

For the sum rules based on f' T fp structures (see Eq.([TT])) we obtain the general results 



F?{t) = -FY{t); F 2 A (t)=F 2 v (t), 



v i 



(77) 



and only the perturbative diagram contributes to F 2 V ^ '. Since the form factors F^^ 
contribute to the semileptonic decays at 0{m 2 /q 2 ) [R| and are thus difficult to measure, we 
will not present results for these quantities here. 
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The stability behavior of the sum rules as a function of the Borel mass is here very similar 
to that shown in Figs. 7,9 and 12, and the pole fits also describe very well the t-dependence 
of the form factors (as well as the curves shown in Figs. 8,10 and 13). In Table IV we give 
the pole parametrization of the form factors for the process A& — > A c lv(_ for two values of k 
(see Eq. Q2"2])) and for the two forms of regions of continuum contribution considered. The 
other parameters are exactly the same as discussed in section IV. The continuum thresholds 
are obtained using method 1 for and A c , and method 2 for A. 



Table IV 

Pole parametrization of the form factors for the process — » A c £f£ for different values of the 

parameters and choices of continuum model. 



continuum model 


K 


FY it) 


F^tXGeV- 1 ) 


rectangular 


1 


7.73/(17.79 -t) 


-0.20/(13.51 - 


rectangular 


2 


8.93/(18.70 -t) 


-0.22/(14.91 - 


triangular 


1 


16.60/(24.24 -t) 


-0.55/(18.63 - 


triangular 


2 


19.02/(25.64 -t) 


-0.61/(19.10- 



We can see from Eq. ([32]) that the IW function at 0{l/m) is given by = —Ff — 
(Ma f — Ma 7 )-F 2 A - The extrapolation of the fits given in Table IV to the maximal momentum 
transfer value yields G^(t maiX ) = 1.0 (0.93) for k = 2 (1) for the rectangular region, and 
G^(t max ) = 1-04 (1.03) for k = 2 (1) for the triangular region, in very good agreement with 
the HQET prediction. 

With the pole parametrizations given above, the A^ — > A c £u£ decay rate is 

T(A b -> Kjrvi) = (5.1 ± 1.8) x 10" 14 GeV , (78) 

where the errors reflect variations of k from 1 to 2 and the different choices for the continuum 
contribution. This value is larger than the one obtained based on // and fp structures 

38 



(see Eq. ([45])), but it is still in agreement with other predictions |16|j27| -p9l and with the 
experimental upper limit [RO] (see Eq. ([46])). 



In the case of the process A c — > A£u£, when the form factors are extracted from the sum 
rules based on f'jfp structure, we obtain a very large decay rate 

T(A+ — > A + e + + z/ e ) = (19 ± 5) x KT 14 GeV , (79) 

which is much larger than the experimental value |3(J (see Eq. (f46"l)) and more than two times 
larger than the value obtained using the sum rules based on fifp structure (see Eq. ([49])). 
The mathematical reason for this discrepancy is that while the F^ A ^ form factors extracted 
from both fifp and f'jfp structures are in very good agreement with each other, this is not 
the case for the form factors, which are zero when the fifp structure is used and of 

about — 0.1(— 0.2) (for k = 1(2)) at t = for the f'jfp structure. This can be seen in Table 
V, where we give the pole parametrization of the form factors for the process A c — > Klvi 
extracted from the sum rules based on the f'jfp structure. 
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Table V 

Pole parametrization of the form factors for the process A c — > htvg for different values of the 
parameters and choices of the continuum model. 



continuum model 


K 


F?(t) 


F 2 v (t)(GeV' 




rectangular 


1 


1.46/(2.98 


-*) 


-0.19/(1.80- 


~t) 


rectangular 


2 


1.78/(3.52 


-*) 


-0.21/(1.94- 


~t) 


triangular 


1 


2.80/(4.10 


-*) 


-0.55/(2.84- 


~t) 


triangular 


2 


3.06/(4.59 


-*) 


-0.53/(3.00- 





The contribution of F 2 is responsible for the strong differences in the results for the 
calculated decay rates. In the A& — ► A c £u£ case, an inspection of Tables IV and V shows 
that although the F^ A ^ form factors are of the same order of magnitude in the A b and A c 
decays, the F 2 ^ form factors are one order of magnitude larger in the A c decay than in 
the A b decay. The results for the A b decay did not change much because the values of the 
form factors F 2 V are consistent with zero for both /j/p and f'jfp structures. 

Appendix 3 . Form Factors in the fif' F Structure 

The observed discrepancy in the results for F 2 V ^ form factors comparing the fifF and 
f'jfp structures for the A c — > ktvi process described above leads us now to consider the form 
factors obtained through the fif' F structure. For these sum rules we obtain 

*t(t) = -F?(t); F 2 A (t) = -F^(t). (80) 

Comparing Eq. (|80|) with Eq. (]77|) we see that the general results are correct only if F 2 V ^ A ' > = 
0, which is true only for the JiJf structure. Therefore, we can foresee more discrepancies 
here. For the F\ form factor, this structure gives a result about 30% smaller than in the 
previous studied structures for the A^ decay, but gives a result with a factor more than two 
smaller for the A c decay. Moreover, while it is still possible to parametrize the t-dependence 
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of Fi with a monopole form in the case of the A& decay (with an agreement similar to the 
other structures), this is not possible for the A c decay. The pole mass needed to fit the t 
behavior is smaller than t max = (M\ c — Ma) 2 . We have used an exponential form factor of 
the type exp (t — a)/b to fit the t-dependence of Fi, but the agreement is rather poor. 

For the F 2 form factor a monopole fit is possible for both decays, and its magnitude is 
about 20% smaller than for the structure f'jfp for the A& decay (therefore, consistent with 
zero), and a factor more than two smaller for the A c decay. 

In Table VI we give the pole parametrization of the form factors for the two processes 
Aft — > A c £u£ and A c — > Aiue, extracted from the sum rules based on the fif' F structure. 



Table VI 

Parametrization of the form factors for the processes A& — > A. c £vg and A c — > A£i/g for different 
values of the parameters and continuum models. 





A b -> AJu e 


A c — > A£z/£ 


continuum model 


K 


F?(t) 


i^WCGeV- 1 ) 


F^(t) 


^(GeV 1 ) 


rectangular 
rectangular 
triangular 
triangular 


1 
2 
1 
2 


6.66/(20.27 -t) 
8.13/(22.50 -t) 
13.74/(26.68 - 1) 
16.17/(29.12 -t) 


—0.21/(15.15 — t) 
-0.22/(13.63 -t) 
-0.41/(18.65 -t) 
-0.45/(19.04 -t) 


e (t-2.07)/0.95 
e (t-2.08)/0.97 
e (t-2.03)/1.09 
e (t-1.72)/1.02 


-0.15/(2.85 -t) 
-0.17/(3.38 -t) 
-0.18/(2.79 -t) 
-0.20/(3.05 -t) 



With the above parametrizations we obtain for the the decay rates 

r(A 6 -> AJ-Pe) = (1.8 ± 0.8) x 10~ 14 GeV, (81) 

and 

T(A+ -> A + e + + z/ e ) = (2 ± 1) x 10" 14 GeV. (82) 

Therefore, the sum rules based on the structure fif' F give too small decay rates for both A& 
and A c decays. Of course the variation in the three structures considered is much larger in 
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the case of the A c decay and this should be expected since the simple current used is not the 
best suited to describe A, as can be seen by the short-dashed line in Fig. 11, which shows 
that the 1 structure is not trustworth in the mass sum rule, and this structure is related 
with the /' structure in the decay. 

It is important to mention that even with this structure the extrapolation of the fits 
given in Table VI to the maximal momentum transfer in the A& decay yields for the IW 
function the value Gf (t max ) = 0.97 ± 0.04 (where Gf = -Ff - (M Ap - M Al )Ff), in very 
good agreement with the HQET prediction. 
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